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>• I Abstract. A class of Laplace transforms is examined to show that par- 

ticular cases of this class are associated with production-destruction and 



(N 



'sj- ■ react ion- diffusion problems in physics, study of differences of independently 

^ . distributed random variables and the concept of Laplacianness in statistics, 

\Q '. a-Laplace and Mittag-Leffier stochastic processes, the concepts of infinite 

O I divisibility and geometric infinite divisibility problems in probability theory 

and certain fractional integrals and fractional derivatives. A number of ap- 
c^ ■ plications are pointed out with special reference to solutions of fractional 

reaction and reaction-diffusion equations and their generalizations. 

%j : 1 Introduction 

We start with simple examples and then will proceed into generalizations. 
Consider a gamma density 

fi{x) = ^^p7^, 2: > 0, /? > 0, Re{a) > 0, 

where Re{-) denotes the real part of (■). Its Laplace transform, with param- 
eter p, is given by the following: 



LfM = / e"PVi(a:)dx = (1 + Pp)-^, 1 + Pp>0. (1) 

Jo 

For a = 1 we have the Laplace transform of the exponential density. Consider 
the Laplace transform of the Laplace density 

1 _M 
/2(a^) = ^e /s , -cx) < X < oo, (3 > 0. 

Its Laplace transform is given by the following 

Lf^ip) = {1 - pvr\ (2) 

This Laplace transform Lf^{p) can also arise from a production-destruction 
or input-output model of the following type: Consider the residual effect of 
an input-output type of situation u = Xi — X2 where Xi and X2 are inde- 
pendently and identically distributed gamma type input and gamma type 
output variables respectively. Then the Laplace transform of the density of 
M, denoted by Lu{p), is given by 



K{p) = L,ML.,{-P) = (1 + /?P)-"(1 - /5p)-" = (1 - /?V)-". (3) 

This residual variable u has applications in many different areas. A few of the 
applications are mentioned in Mathai (1993a) and in Mathai, Provost, and 
Hayakawa (1995). The Laplace transform in (3) is also associated with the 
concept of Laplacianness in statistics, see Mathai (1993). This concept is con- 
nected to bilinear forms, quadratic form and the concept of chi-squaredness 
of quadratic forms, which is the basis for making inference in analysis of 
variance, analysis of covariance, regression and general model building areas. 
If there are several independent input variables xi,---,Xn such as the 
situation in reaction or production problems, and if there are several inde- 
pendent output variables Xm+i,- ' ' ^^m+n and if they are all gamma type 
variables with different parameters then the residual u = xi + ■ ■ ■ + Xm — 
Xrn+i — ■ ■ ■ — Xm+n has the Laplacc transform 



L„(p) = (l + /5ip)--i...(l + /?„p)-» 

X (l-/3„+ip)-"-+^---(l-/9„+„p)-"-+". (4) 



The density of u, corresponding to the Laplace transform in (4), is also 
available in the literature, see for example Mathai and Provost (1992). Now 
consider the Mittag-Leffler function 
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^^r{l + ka) 
and let 

/3(x)=i?„(-x"), 

then the Laplace transform 

Lf.^{p)= e-P^£;„(-x")dx=p"-i(l+p")-\ (5) 

JO 

Haubold and Mathai (1995, 2000) considered a reaction problem where 
the number density of the reacting particles is a function of time t. For the 
i*^ particle let the number density be Ni = Ni{t) with Ni{t = 0) = A'^o- If 
the production rate is proportional to the number density then we have the 
simple differential equation, dropping i, 

AiV(t) = a Nit), 

where a is a constant. If some particles produced are also destroyed or 
consumed and if the destruction rate is given by 

AiV(t) = -6 N{t), 

then the residual effect is given by 

d 

—N(t) = -c N(t), c=b-a, 
at 

and then the solution is 

N{t) = A^oe"'*. 

But if a fractional integral, instead of the full integral, is used then the 
reaction equation is given by 

N{t) -No = -c\Dr^N{t),v > 0, (6) 

where c is replaced by c'^ and its solution has the form 



where the fractional integral is defined as follows 

aDr-'fit) = -^ fit - ur-'f{u)du, z. > 0, (7) 

with aDt* f{t) = f{t). The Laplace transform of N(t), coming from (6), is 
then, 



AT, 
p[l 

which is a constant multiple of (5) for c = 1,1/ = a. 



LNit)ip) = zjrT^^^ (8) 



2 Generalizations 

If the basic fractional production-destruction equation in (6) is modified to 
the form 

N{t) - Nof"-^ = -c\Dt-''N{t),ii > 0, z/ > 0, (9) 

then the Laplace transform is given by 

tor II = V. (llj 



This Laplace transform in (11) for c = 1 is associated with infinitely divisible 
and geometrically infinitely divisible distributions, a-Laplace and Linnik dis- 
tributions in probability theory and a-Laplace and Mittag-Leffler stochastic 
processes, see the works of R.N. Pillai and his associates, a summary of which 
is available from Jose and Seetha Lekshmi (2004). A comprehensive presen- 
tation of this field of processes is given in Uchaikin and Zolotarev (1999). 
The class of Laplace transforms relevant in geometrically infinitely divisible 
and a-Laplace distributions is of the form 



where r]{p) satisfies the condition 'r]{bp) = h'^r]{p) where rj{p) is a periodic 
function for fixed a. Observe that (11) satisfies the above condition with 
rj{p) = p^ and a = v. Thus the basic reaction equation in (9) is also associ- 
ated with certain stochastic processes and infinitely divisible distributions. 

Another observation that can be made is that the inverse of (10) or the 
number density N{t), which gives rise to (10), is of the following general form 

N{t) = N.Vi^Y-^E.^^i-d'n, (13) 

where Ev^^{-) is a generalized Mittag-Leffler function, 

oo ( 1 \kfj.v „v\k 

^..(-«T-g^^(^./'>o.">o. (14) 

Note that (14) is also connected to the F-function introduced by Hartley 
and Lorenzo (1998) and Lorenzo and Hartley (1999). Their F-function is the 
following, 

F,{-a,t) = f-' Y: L i ^ = t'-'E,,,i-an,Reiq) > 0. (15) 
^oT{q + nq) 

When the variable is restricted to the interval t > 6 for some S, then (15) 
reduces to the R-function of Lorenzo and Hartley (1999), which is defined as 



R,,^{a,S,t) = J2 r\/ Jn ^' t>S>0 

^0 T[{n+l)iy-fi] 

= {t-dy-^'-^E.^.^^lait-SY], t>6>0. (16) 

If we consider another modification of the basic reaction equation in (9) to 
the form 

N{t) - No f^-'El^i-c^n = -c\Df-N{t), (17) 

where E^ J-) is a further generalized form of the Mittag-Leffler function given 
by 



(7)/.- = 7(7 + 1) • ■ ■ (7 + fc - 1), 7 # 0, (7)0 = 1, 
then the Laplace transform is the following 

The inverse of (19) or the number density in this case is given by 

N{t) = No t^'-^El+^^-d't^). (20) 

The most general form of Laplace transform associated with Mittag-Leffler 
functions is that in (19), which is a special case of the general class of Laplace 
transforms associated with a-Laplace stochastic processes and geometrically 
infinitely divisible statistical distributions. 

If the Mittag-Leffler function is further generalized then the Laplace trans- 
forms will enter the family of Fox's H-function. For example, let us consider 
a simple generalization of the Mittag-Leffler function to the following form 

^^ t'o(f^^)^k\ T{f3+ka)- 

Then the Laplace transform of t^~^ giit'^) , denoted by Li{p), is given by the 
following 

LiiP) = E 7^ P^^"'^^^ =p-^Fi(7i; A;p-"), (21) 

where iFi is a confluent hypergeometric function. If the factor t^~^ is not 
incorporated then the Laplace transform will be a special case of an H- 
function. For a description and properties of Fox's H-function see Mathai 
and Saxena (1978), Mathai (1993a), and Kilbas and Saigo (2004). 

A direct generalization of a Mittag-Leffler function and the most gener- 
alized form in this category is Wright's function given by 



i^q 



pYq 



r(ai,Ai),...,{ap,Aj,) n ^ {l\%iT{aj + Ajk)} z'' 

' z — 
{h,Bi),-,(bg,Bg) ' . 



to{uUnbj + B,k)}k\' 



with 1 + E'=i B, - Ef=i A, > 0. 



Wright's function is again a particular case of Fox's H-function. Note that 
the Mittag-Leffler function 



Ea,/3iz) = 1^1 



r{i.i) 



; z 






I (0,1) 
l(0,l),(l-/3,«) 



(23) 



where -f/^(-) denotes an H-function, and the H-function is defined in terms of 
the following Mellin-Barnes type representation 



TTm,n 



<{a-i,Ai),--,(ap,Ap) 
\(bi,Bi),-,ibq,Bq) 



1 

2ni Jl 



g{s)z ^ds,i 



-1, 



where 



9{s) 



{nr=i r{b, + B,s)} {uu r(i - % - A,s)} 

{nj=™+i r(l - b, - B,s)} {U%n+i r(a, + A,s)} 



(24) 



(25) 



with Aj,j = 1,- ■ ■ ,p and Bj,j = 1,- ■ ■ ,q being positive real numbers and o/s 
and fej's are complex quantities. Existence conditions, the types of contour 
L and properties are available from Mathai and Saxena (1978) and Kilbas 
and Saigo (2004). 

Observe that for 7 = 1 in (18) we have El^^ = E^^^. Hence for 7 = 1, 2, 3 ■ ■ ■, 
for positive integer values, one can write i?J as a linear function of £"., .(■). 
Thus, certain linear functions of Laplace transforms of the type in (10) can 
be summed up to a Laplace transform of the type in (19). Some examples 
of this type may be seen from a series of papers by Haubold, Mathai, and 
Saxena, for example, see Saxena, Mathai, and Haubold (2004). One such 
example can be stated as follows: Ifc>0,yU>0,z/>0 then the solution of 
the fractional integral equation for the number density N{t), 

Nit) - No t^-'E.^^i-c-'n = -c\Dt~''N{t), (26) 



E.^^^.i-c^n + (1 - ^ + z/)E,,^(-c'^r) . (27) 



is given by 




Nit)- 


_ No t^-i 


V 



3 Reaction Equation 

If the Laplace transform has the structure of a product of various factors 
then such a case can also be handled without much difficulty. As an example, 
consider the production-destruction fractional integral model for the number 
density N(t) of the following type 

N{t) - No t'^-'E^^^i-d-'r) = -c\Df-N{t). (28) 

Then it is not difficult to see that the Laplace transform of N{t), with pa- 
rameter p, denoted by L2{p) is the following 



L2{p) 



No 



(29) 



P^[1 + (P1[1 + (P1- 

When c = d then this case reduces to the Laplace transform of a generalized 
Mittag-Leffler function, which can also be written as a linear combination of 
simple Mittag-Leffier functions as explained in (26) and (27). When c ^ d 
then we may consider the following identity 



{p" + C) {p" + d") {c" - d") 



p"" + d" p" + & 



,c^d. 



Hence, 



(30) 



L2{p) 



1 



1 



1 



1 



1 



{d" - d") I p^'~'' 1 + {^Y pf-'' 1 + {^y 



Then taking the inverse Laplace transform, we have. 



N{t) = Ni 







c 



d" 



E i 



-d^e 



E ( 



-cH" 



(31) 



(32) 



Now, let us look into another situation where the Laplace transforms are 
of the following types 



LM 



p 



n-l 



and Li{p) 



p 



,/3-l 



(33) 



where a, 6, a, (3 are constants. Such Laplace transforms can be handled by 
using the following procedure. For example, consider 



Lsip) 



„o-l „a-/3-l 



oo -(a-/3)r-l 



r=0 



Comparing with the Laplace transform of a generahzed Mittag-Leffler func- 
tion we have the inverse given by the following 



L-^ 



^a-l 



p" + apl^ + h 



Y.i-^n^''-^^''Kt^p)rA'bn- (35) 

r=0 



4 Reaction-Diffusion Equation 

A problem recently considered by Saxena, Mathai, and Haubold (2005) is a 
react ion- diffusion system occurring in various areas of pattern formation in 
biology, chemistry, physics, see for example Henry and Wearne (2000, 2002), 
Henry, Langlands, and Wearne (2005), and Manne, Hurd, and Kenkre (2000). 
The basic reaction-diffusion equation has the form 

dN d'^N 

^ = f.^ + \f{N),N = N{x,t), (36) 

where /i is the diffusion constant, A is a constant and f{N) is a nonlinear 
function of A^. When f{N) = 6N{1 — N), where 5 is a constant, then (36) 
reduces to the real Ginsburg-Landau equation. An equation in this category 
examined by Manne, Hurd, and Kenkre (2000) is the following 

where ^ indicates the strength of nonlinearity in the system. For solving (37) 
one can adopt the procedure of taking the Laplace transform with respect to 
t and Fourier transform with respect to x, simplifying and then inverting the 
Laplace- Fourier transform to obtain N[x,t). Details of the procedure may 
be seen from Saxena, Mathai, and Haubold (2004, 2005). In this procedure, 
the crucial step is the inversion of the Laplace transform of the type L^lp) 



in (33). From (35) one can see that the inverse can be written as a series of 
Mittag-Leffler functions 



p/3-l 



(_«)^,(-/^)(^+i)£;^+i (_5,«). (38) 



Z^V "; f- ^a,(a-/3)(r+l)+l 



.p" + api^ + b. 

5 Conclusions 

Linear and nonhnear reaction, diffusion, and reaction-diffusion equations, 
respectively, are used to model spatio-temporal processes in physical sys- 
tems, for example astrophysical fusion plasmas (Wilhelmsson and Lazzaro, 
2001; Kulsrud, 2005) . Such equations give rise to dissipative structures 
and self-organization phenomena (Nicolis and Prigogine, 1977; Haken 2004). 
Attempts have been made to employ them to discover a time-dependent 
mechanism in solar nucleosynthesis (Haubold and Mathai, 1995, 2000). 

The linear reaction (relaxation) equation can be used to explore funda- 
mental principles of Boltzmann-Gibbs statistical mechanics and its nonlinear 
generalization leads to new insights into nonextensive statistical mechanics 
(Tsallis, 2004). This approach has also been extended to discover the con- 
nection between Tsallis' nonextensive maximum entropy formalism and the 
nonlinear reaction-diffusion equation (Tsallis and Bukman, 1996). 

In this paper we used Laplace transform techniques to derive closed-form 
solutions of reaction equations starting from its simplest form and its frac- 
tional version in eq. (6) through generalizations of the equation given in 
(9), (17), (26), and (28). Along this way of generalizing the linear reaction 
equation we provide their respective closed-form representations in terms 
of Mittag-Leffler functions. Gradually the link between the Mittag-Leffler 
function and Wright's function and, subsequently. Fox's function has been 
shown. This transition can be characterized as a transition from the stan- 
dard exponential behavior ("normal" relaxation) to a power-law behavior 
("anomalous" relaxation) of the solutions of the respective reaction equa- 
tions. For the react ion- diffusion equations (36) and (37) we indicate that 
the same Laplace transform techniques, as successfully employed for deriv- 
ing closed-form solutions of reaction equations, can also be used for finding 
solutions of reaction-diffusion equations. 
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